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J/ip production in d + Au and p + p collisions at ^snn = 200 GeV has been measured by the 
PHENIX experiment at rapidities —2.2 < y < +2.4. The cross sections and nuclear dependence of 
J/tp production versus rapidity, transverse momentum, and centrality are obtained and compared 
to lower energy p + A results and to theoretical models. The observed nuclear dependence in d + Au 
collisions is found to be modest, suggesting that the absorption in the final state is weak and 
the shadowing of the gluon distributions is small and consistent with Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi-based parameterizations that fit deep-inelastic scattering and Drell-Yan data at lower 
energies. 
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J/tp production in hadron collisions, since it proceeds 
predominantly through diagrams involving gluons (e.g. 
gluon fusion) |l(, is a sensitive probe of the gluon struc- 
ture function in the nucleon and its modification in nu- 
clei. It is also a leading signal for the creation of hot- 
dense matter in heavy- ion collisions |2|. Shadowing of 
partons (quarks or gluons) in nuclei is a depletion of 
their population at small momentum fraction of the nu- 
cleon, x, compared to that in a free nucleon, with a corre- 
sponding enhancement at moderate x (anti-shadowing). 
In y/sNN = 200 GeV deuteron-gold (d + Au) collisions at 
the Relativistic Heavy Ion Collider (RHIC), for positive 
(deuteron direction) rapidities, gluons are probed that lie 
well into the shadowing region with momentum fractions 
in Au, X2 ~ 3 x 10~ 3 . Models of gluon shadowing predict 
suppressions of J/ip production in nuclei that differ by as 
much as a factor of three 0, [3- • Recent theoretical de- 
velopments, e.g. the Color Glass Condensate model 0, 
suggest that at very low x non-linear gluon saturation 
effects become important and cause substantial modifi- 
cations of the gluon densities. 

The connection of the observed J/tp suppression to the 
modified gluon distribution in nuclei can be clouded by 
the absorption of the final-state cc which depends on 
the poorly known production mechanism and by the 
energy loss of the intial-state gluon - although the latter 
is thought to be small at RHIC energies yj. This con- 
nection is also distorted by the fact that approximately 
a third of the J/tp 's come from decays of higher-mass 
resonances yf . 

Here we present measurements made by the PHENIX 
experiment at RHIC for the production of J/tp 's in 
tJsnn — 200 GeV d + Au and proton-proton (p + p) 
collisions. These data provide the first measurement of 
the nuclear dependence of J/tp production at this en- 
ergy, a much higher energy than previous p + A mea- 
surements from fixed-target experiments at ^/snn 40 
Although our measurements are 
for d + A, the nuclear effects on the J/tp in deuterium 
were found to be small at lower energies yT3|. Besides 
the shadowing region at small x, these data also probe 
larger gluon momentum fractions (at negative rapidity) 
nearer the rest frame of the residual nucleus. Finally, 
these measurements also serve as a baseline for the up- 
coming results from the high-luminosity Au + Au and 
Cu + Cu runs and must be understood in order to look 
for effects beyond what is expected from cold nuclear 
matter. 

The measurements described here are similar to earlier 
ones with PHENIX for p + p HH and Au + Au H3 
collisions, but with a second muon spectrometer added 
and higher luminosity. The two muon spectrometers 
are especially valuable for asymmetric collisions such as 
d + Au where simultaneous measurements at positive 
(1.2 < y < 2.4) and negative (-2.2 < y < -1.2) rapidi- 
ties, along with central (\y\ < 0.35) rapidity from e + e~, 



are then available. Electrons in the central arms are iden- 
tified by matching charged particle tracks to clusters in 
an electromagnetic calorimeter (EMC) and to rings in 
a ring imaging Cerenkov (RICH) detector. Muons are 
identified by their detection in Iarocci tubes after their 
penetration through 8 to 11 interaction lengths of copper 
and steel absorber. 

The data used in this analysis were recorded in 2003 
using a trigger that required hits in each of the two beam- 
beam counters located at negative and positive rapidity 
(3 < \ri\ < 3.9). In addition, for the di-muons at least 
two tracks in the muon identifier of appropriate absorber 
depth were required, while for the di-electrons a one-track 
trigger with a signal above threshold in the EMC with 
a matching hit in the RICH was required. After quality 
and vertex cuts, the samples for the three arms corre- 
spond to integrated luminosities from 180 to 250 nb _1 
(p + p) and 1.4 to 1.7 nb _1 (d + Au), 

For the di-muons the J/tp yield is obtained after sub- 
traction of the combinatoric background using like-sign 

muon pairs (2^J N ++ N ) and by fitting the resulting 

mass peak with a Gaussian plus an exponential to rep- 
resent the small remaining continuum background un- 
derneath the peak. A variety of continuum shapes were 
checked for each fit in order to establish the uncertainty 
due to the low-statistics background. For the di-electrons 
the combinatoric background was subtracted using the 
sum of like-sign pairs and the J/tp yield was taken as all 
remaining events in the mass range 2.6 to 3.6 GeV/c 2 . 
A total {p + p plus d + Au) of about 2100 and 500 J ftp's 
were obtained in the /i/i and ee channels, respectively. 

The differential cross sections are calculated as, 
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and the nuclear modification factor, RdA, is 



RdA = 



dajfpdy 



(2 x 197) x daf^/dy 
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In the above expressions Bu is the J/tp branching ratio 
to di-leptons, Nj/$ is the measured J/tp yield, A is the 
geometrical acceptance, e rec is the di-lepton reconstruc- 
tion efficiency, etrig is the trigger efficiency, N evt is the 
number of min-bias triggers sampled, o M l B is the total 
minimum-bias (MB) cross section, Ay is the rapidity bin 
width, and e M B B c and Ejffi are the beam-beam trigger 
efficiencies for min-bias and J/tp events respectively. The 
factor of 2 x 197 causes RdA to be one if the d + A cross 
section is just additive from p + p, i.e. if there are no 
nuclear modifications. 

For p + p we use the cross section for our beam-beam 
trigger, cr% t B {pp)e M B B c {pp) = 23.0 ± 2.2 mb; and the effi- 
ciency for events with a J/tp, £ c ' (pp) = 0.79±0.02. For 
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d+Au collisions we use a beam-beam trigger cross section 
(dAu) = 1.99 ± 0.10 b from our mea- 



of alf B (dA£ 



C BBC 
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surement [1£| using photo-dissociation of the deuteron 
as a reference which is consistent with our calcu- 
lated Glauber result of 1.92 ± 0.18 b. For the J/tp we use 



r BBC 

-J/i> 



(dAu) = 0.94 ±0.02 [2C 



For d + Au collisions, the centrality of the collision can 
be characterized by measuring the charge deposited in 
the beam-beam counter in the Au beam direction [20| . 
An approximate number of nucleon+nucleon collisions 
(N co u) can be obtained through a Glauber calculation 
that relates this (N co u) to the observed charge. In this 
case RdA{N co ii) is calculated as, 
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where the invariant yield N inv is, 



N inv ((N coH )) 



Nj/^Cbias((N C oll)) 
A tree ttrig [N ev t ( Aw / w)] 



(4) 



with (N co ii) being the average number of binary collisions 
for a particular d + Au centrality bin and N evt (Aw/w) 
the number of d + Au min-bias triggers sampled that lie 
in this fraction, Aw, of the total minimum-bias centrality 
range, w. This prescription is equivalent to that of Eq. 
1,2 for minimum bias. For p + p collisions Aw/w is one. 
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FIG. 1: (color online) (a) The 200 GeV J/tp p + p differential 
cross section times di-lepton branching ratio versus rapidity 
(10% overall normalization uncertainty is not included), (b) 
The minimum bias RdA versus rapidity (12% overall normal- 
ization uncertainty not included). For both panels the dashed 
error bars represent systematic uncertainties relevant for com- 
paring the two rapidity bins in each muon arm, while the solid 
error bars represent the overall uncertainties relevant for com- 
paring points at negative, central, or positive rapidity. The 
curve in (a) represents a fit as described in the text while 
the curves in (b) are theoretical calculations 0, l24l |2^| as 
described in the text. 
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/(.bbc j g a correc tion for the smaller trig- branching ratio of 5.9% 23], gives a total cross section 



ger efficiency in minimum-bias events compared to those 
with a J/ip- For d + Au, Cuas depends on (N co u) and 
takes into account the effect of the underlying event mul- 
tiplicity on both the trigger efficiency and the centrality 
measurement |20j . Its variation with centrality is up to 
7% from unity. 

For the electron analysis, Ae rec and t trig were deter- 
mined using a Geant [21| simulation of the central arms 
and a trigger response software emulation 0|. e rec was 
confirmed by studying pairs identified as photon conver- 
sions in the data. The systematic uncertainty of 10.4% 
is dominated by run-to-run efficiencies (5%), yield ex- 
traction (5%), and the occupancy dependence of the ef- 
ficiency (4.4%). 

For the muon arms, A e rec ttrig was determined within 
each rapidity and pt bin, using a Geant simulation with 
J/ip events generated by Pythia . The dominant sys- 
tematic uncertainties in our result arc +6/— 9% from the 
muon identifier efficiency and up to 10% (for the most 
central negative rapidity d + Au data) from the combi- 
natoric background. 

Figure^(a) shows the measured p+p differential cross 
section times branching ratio versus rapidity with a di- 
electron point at mid-rapidity and two di-muon points 
at negative and positive rapidities. A fit to a shape gen- 
erated with Pythia is performed and, using a di-lepton 



j-p^ = 2.61 ± 0.20(fit) ± 0.26(abs) fib. Variations in the 
parton distribution functions and models used to deter- 
mine the shape are negligible compared to the fit errors. 
This result is smaller by about two sigma than our pre- 
vious lower statistics result fl6| . 

Figure^] (b) shows the nuclear modification factor RdA 
(Eq. 2) versus rapidity, where a value of one would cor- 
respond to no nuclear modification. While this ratio is 
consistent with unity at negative rapidity, it is signif- 
icantly lower at the most positive rapidity where glu- 
ons are expected to be shadowed in a heavy nucleus. 
Theoretical predictions 0, 0> El] that include the ef- 
fects of absorption and shadowing are shown for com- 
parison in Fig. H (b) . The data favor a relatively mod- 
est shadowing in agreement with the parametrization 
of Eskola-Kolhinen-Salgado (EKS) based on a leading- 
twist Dokshitzer-Gribov-Lipatov-Altarelli-Parisi-evolved 
parametrization of nuclear deep inelastic scattering and 
Drell-Yan data at lower energies rather than the 
stronger gluon shadowing of Kopeliovich || or Frankfurt- 
Guzey-Strikman (FGS) [4| based on models involving co- 
herence for a qq dipole in the nucleus. The cc absorption 
cross section is not well determined by our data, but is 
probably nearer to 1 mb and is certainly smaller than the 
4.1 ± 0.4 mb found at lower energy |26j. 

Lower energy p + A measurements showed that J 
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FIG. 2: (color online) a versus (a) a versus pp compared 
to lower energy measurements shown for three different xf 
ranges. An additional overall uncertainty of 0.02 in our a 
values is not shown. 



suppression did not follow a universal behavior vs X2 |27fl , 
the momentum fraction in the heavy nucleus, as would 
be expected if the suppression was dominated by shad- 
owing [icj . As shown in Fig. (a), our data confirm 
this X2 scaling violation with the addition of a smaller x-i 
point, but with our small range in a;;? have little to add 
to the approximate xp scaling observed in these lower 
energy measurements [see Fig. 121(b)]. Here a is defined 
by <JdA — &pp x (2^4) Q and xp = x\ — X2, where X\ is the 
momentum fraction of the gluon in the deuteron. This 
xp scaling may be caused by energy loss of the gluon 
in the intial state 0, or by an energy conservation ef- 
fect ( "Sudakov suppression" ) which causes a univer- 
sal suppression that increases with xp for production of 
J/ip's and other hadrons. It is also similar to that ob- 
served for more positive rapidity hadron production in 
the d + Au "limiting fragmentation" region [2^| . 

Invariant cross sections versus transverse momentum, 
(cP a I dydpr ) /(2irp T ) , have been fit to the form A x (I + 
{pT I B) 2 )~ 6 [3fJ. Average p T values resulting from these 
fits are 4.28 ± 0.31, 3.03 ± 0.40 and 3.63 ± 0.25 GeV/c 2 
for d + Au collisions at negative, zero, and positive xf, 
respectively; compared with 2.51 ± 0.21 and 4.31 ± 0.85 
GeV /c 2 for negative/positive and zero xp p+p collisions, 
respectively. The observed pp broadening is shown in 
Fig. (a) . For negative xf it is consistent with that of 
the lower energy (^Jsnn = 39 GeV) measurements from 
E866/NuSea |l0|. but may be flatter at positive xf- At 
zero xf no pp broadening is seen within errors. 

RdA (Eq. 3) is shown in Fig. |3] (b) for four centrality 
classes and for minimum bias collisions. This classifi- 
cation into centrality bins for these results can only be 
approximate, as indicated by the overlapping histograms 
of N co u. At positive rapidity (small X2, or the shadowing 



region), a weak drop for more central collisions is ob- 
served, while no significant centrality dependence is seen 
for negative rapidity or for central rapidity. The theoreti- 
cal curves on Fig. [21(b) correspond to different amounts of 
density dependent shadowing and anti-shadowing 0,113 
and also include absorption. They are consistent with our 
data except at positive rapidity where the EKS shadow- 
ing curve is closest to our results, although slightly lower 
perhaps due to the amount of absorption that is included. 

In summary, during the RHIC 2003 run, the PHENIX 
experiment measured nuclear effects on J/-0 production 
for d + Au collisions at ^/snn — 200 GeV. Increasing 
suppression for larger rapidity (smaller X2) and for more 
central collisions (higher nuclear densities sampled) both 
are consistent with models containing a small amount of 
impact-parameter dependent shadowing and with weak 
absorption. Theoretical calculations which include EKS 
shadowing seem most consistent with the data. However 
comparisons with other measurements at lower energies 
show that shadowing cannot be the dominant effect, at 
least not for the lower energy measurements. We also see 
some transverse momentum broadening which is consis- 
tent with that seen at lower energy. Higher luminosity 
d + Au running in the future yielding higher numbers of 
J/tjj 's will be necessary to quantify these nuclear effects 
and to more clearly distinguish between various theo- 
retical models of shadowing, absorption and other cold 
nuclear matter effects. 

We thank the staff of the Collider-Accelerator and 
Physics Departments at BNL for their vital contri- 
butions. We acknowledge support from the Depart- 
ment of Energy and NSF (U.S.A.), MEXT and JSPS 
(Japan), CNPq and FAPESP (Brazil), NSFC (China), 
IN2P3/CNRS, CEA, and ARMINES (France), BMBF, 
DAAD, and AvH (Germany), OTKA (Hungary), DAE 
and DST (India), ISF (Israel), KRF and CHEP (Ko- 
rea), RMIST, RAS, and RMAE (Russia), VR and KAW 
(Sweden), U.S. CRDF for the FSU, US-Hungarian NSF- 
OTKA-MTA, and US-Israel BSF. 



* Deceased 

' PHENIX Spokesperson zajc@nevis.columbia.edu 
[1] M. Bedjidian et al. hep-ph/0311048 (2003). 
[2] T. Matsui and H. Satz, Phys. Lett. B178, 416 (1986). 
[3] K. J. Eskola, V. J. Kolhinen, and R. Vogt, Nucl. Phys. 

A696, 729 (2001). 
[4] L. Frankfurt, V. Guzey, and M. Strikman, Eur. Phys. J. 

A5, 293 (1999). 
[5] B. Kopeliovich, A. Tarasov, and J. Hufner, Nucl. Phys. 

A696, 669 (2001). 
[6] L. McLerran and R. Venugopalan, Phys. Rev. D 49, 2233 

(1994). 

[7] R. Vogt, Phys. Rev. C 61, 035203 (2000). 

[8] L. Antoniazzi et al. Phys. Rev. Lett. 70, 383 (1993). 

[9] D. M. Aide et al., Phys. Rev. Lett. 66, 133 (1991). 



6 




FIG. 3: (color online) (a) a versus pr compared to lower energy measurements, shown for three different xf ranges. The error 
bars have the same meaning as in Fig. An additional 0.02 overall uncertainty is not shown. The dashed curves are simple 
fits EH to the lower energy results, (b) Nuclear modification factor versus centrality as given by the number of nucleon+nucleon 
collisions shown for three different rapidity ranges, compared to theoretical calculations|24l. |2E| including final-state absorption 
and EKS (solid) or FGS (dashed) shadowing. The bars at the low end of each plot represent the systematic errors between 
different rapidity ranges. An additional 12% global error bar is not shown. The histograms at the bottom of the lower panel 
indicate the distribution of the number of collisions for each of the four centrality bins. 



[10] M. J. Leitch et al., Phys. Rev. Lett. 84, 3256 (2000). 
[11] B. Alessandro et al., Phys. Lett. B553, 167 (2003). 
[12] J. Spengler J. Phys. G30, S871 (2004). 
[13] J. Badier et al., Z. Phys. C20, 101 (1983). 
[14] M. J. Leitch, Eur. Phys. J A S19, 129 (2004). 
[15] S. S. Adler et al., Nucl. Instrum. Meth. A499, 469 
(2003). 

[16] S. S. Adler et al., Phys. Rev. Lett. 92, 051802 (2004). 

[17] S. S. Adler et al., Phys. Rev. C 69, 014901 (2004). 

[18] S. N. White, Proceedings of the Xlllth International 

Workshop DIS'05, Madison, WI, 2005, edited by W.H. 

Smith [AIP Conf. Proc, 792, No. 1, 527, (2005)]. 
[19] S. R. Klein and R. Vogt, Phys. Rev. C 67, 047901 (2003). 
[20] S. S. Adler et al., Phys. Rev. Lett. 94, 082302 (2005). 
[21] F. Carminati et al., "GEANT: Detector Description and 



Simulation Tool", CERN Program Library Long Writeup 
No. W5013, 1993 (unpublished). 
[22] T. Sjostrand et al., Comput. Phys. Commun. 135, 238 
(2001). 

[23] S. Eidelman et al., Phys. Lett. B592, 1 (2004). 
[24] S. R. Klein and R. Vogt, Phys. Rev. Lett. 91, 142301 
(2003). 

[25] R. Vogt, Phys. Rev. C 71, 054902 (2005). 
[26] G. Borges, et al. hep-ex/0505065 (2004). 

[27] x 2 = 0.5(-2;f + J^+AM^Jsmn). 

[28] B. Z. Kopeliovich, et al. hep-ph/0501260 (2005). 

[29] B. B. Back et al., Phys. Rev. Lett. 93, 082301 (2004). 

[30] J. K. Yoh et al., Phys. Rev. Lett. 41, 684 (1978). 



